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.1, Abstract 

Scaling and hyperscaling laws provide exact relations among critical exponents describing the 
^ behavior of a system at criticality. For nonequilibrium growth models with a conserved drift there 

exist few of them. One such relation is a + z = 4, found to be inexact in a renormalization group 
^ calculation for several classical models in this field. Herein we focus on the two-dimensional case 

o 

^ and show that it is possible to construct conserved surface growth equations for which the relation 

a + z = 4 is exact in the renormalization group sense. We explain the presence of this scaling law 
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if) in terms of the existence of geometric principles dominating the dynamics. 
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I. INTRODUCTION 



One of the great successes of equilibrium statistical mechanics is the classification of 
critical phenomena into universality classes. These are characterized by sets of critical 
exponents which describe the system long range physics in the critical point. It is precisely 
the infrared behavior, which is independent of the microscopic details of the particular 
system, which lies at the basis of universality. 

The numerical values of the critical exponents can be calculated exactly in a number of 
fortunate cases, otherwise they have to be calculated approximately using some perturba- 
tive technique, like notably the renormalization group, or by means of numerical simulations. 
Scaling and hyperscaling laws provide exact relations among critical exponents and conse- 
quently both are a deep expression of the physics and provide a useful calculational tool. 

Nonequilibrium statistical mechanics has been built in many situations as an extension of 
its successful equilibrium counterpart. The classification in terms of universality classes and 
the tools for computing critical exponents have been adapted to systems out of equilibrium, 
such as absorbing state transitions [Ij and stochastic growth [2] . This last field is sometimes 
regarded as paradigmatic within nonequilibrium statistical mechanics. The reason for this 
is the ubiquity of certain universality classes which were originally discovered in this context 
and were subsequently found in different areas of physics. A particularly relevant example of 
this is the Kardar-Parisi- Zhang (KPZ) equation [3J, which scaling behavior has been related 
to phenomena as diverse as far from equilibrium growth, turbulence and directed polymers. 

The KPZ and other nonequilibrium growth equations can be characterized by three scal- 
ing exponents: the roughness exponent a characterizing the interface morphology, the dy- 
namic exponent z which specifies the velocity at which correlations propagate and the growth 
exponent /3 describing its short time dynamics [2]. One of the most important relations 
among these exponents is the scaling relation a = Pz, which is fulfilled in a large number 
of cases. Together with it there exist other scaling and hyperscaling relations describing 
fundamental characteristics of different physical situations. 

This work is devoted to the study of stochastic partial differential equations with a 
conserved drift. By this we mean that the deterministic part is conservative or, in other 
words, the spatial average of the deterministic part of the equation vanishes provided suitable 
boundary conditions are used. 
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In the following we will examine one of this scaling relations already found in the liter- 
ature, and we will investigate its possible physical origin. In this work we concentrate on 
the two-dimensional situation. The paper is organized as follows: in Section II we present 
some scaling relations for nonequilibrium models, in Section III we give a short geometric 
derivation of a family of conserved growth equations we will be interested in. Section IV is 
devoted to the renormalization group analysis of these and other related equations. In the 
last two Sections we discuss the possible origin of scaling relations for the different models 
and draw our conclusions. 

II. SCALING RELATIONS 

The KPZ equation reads [3] 



where ^ is a zero mean Gaussian noise which correlation is given by {C,{r,t)^{f' ,t')) = 
D5{f — f')6{t — t'). If one performs the dilatation r — > 6r, t b^t and h ^ b"' h the 
nonlinearity becomes scale invariant provided a + z = 2. This scaling relation is actually 
fulfilled by the KPZ equation in any dimension and its exactness is usually attributed to 
its Galilean invariance property [H |5]. However the role of Galilean invariance was put into 
question in a related nonequilibrium model [6]-|8]. In particular, in these works it is shown 
that neither Galilean invariance nor extended Galilean invariance contribute to vertex non- 
renormalization beyond the zeroth mode. Also, the same scaling relation a + z = 2 was 
found in numerical schemes |9HTT] and stochastic equations fi2[ 113] which do not explicitly 
fulfill the Galilean invariance symmetry. 

Together with the KPZ equation, in this field one finds the Sun-Guo-Grant (SGG) [H] 
and Villain-Lai-Das Sarma (VLDS) [151 US] equations, which read 



where the noise ^^^^ corresponds to the SGG equation and ^^'^^ to the VLDS one. 

In both cases the noise is a zero mean stochastic Gaussian process and the respective 
correlations are: 



dth = uV^h+-{Vhy + ar,t), 



(1) 



dth = -uV^h + XV^Vhf + ^'~^'^\f,t), 



(2) 



(e(^)(r,t)e(^)(f',0) = -DV'6{f-f')6{t-t') 



(3) 
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and 



f-r')5{t-t'). 



(4) 



For these equations the nonhnearity becomes scale invariant to the dilatation transformation 
provided the exponents fulfill the scaling relation a + z = A. This relation was found to 
be fulfilled at one loop order renormalization group calculations and it was conjectured to 
be exact due to the existence of a functional Galilean invariance symmetry [14J. However, 
in [TTj, Janssen showed that the mathematics leading to this symmetry was ill-posed. Also, 
by means of a renormalization group calculation, he showed that this relation does not hold 
at two loops [T7]. The correction to the scaling relation he found is however small and 
presumably not detectable in simulations or experiments. 

We also would like to mention here that this correction has not been found in the results 
obtained by means of other methods, such that the self-consistent expansion [T8H20], when 
applied to these equations |21J. This approximation scheme has proved itself very useful 
in finding exponent values in some cases in which the renormalization group fails, like for 
instance for the nonlocal counterpart of some stochastic growth equations |22l |23] . On the 
other hand, renormalization group analyses have been able to yield exact results for some 
stochastic growth equations [211 |2S]. We will restrict ourselves to the renormalization group 
analysis of local equations and leave different approaches for a more comprehensive study. 
As will be evident in the following, this analysis has a clear interpretation in the cases under 
consideration. 

In the next section we will show that it is possible to find models for which the scaling 
relation a + 2; = 4 is exact to all orders in a renormalization group calculation. 

III. GEOMETRIC DERIVATION OF A STOCHASTIC GROWTH EQUATION 

We now briefiy review the geometric derivation of a stochastic growth equation [26]. We 
concentrate on d = 2 and assume the following variational approach for the dynamics 



where we have added the noise term .^^^'^^ and X is a nonequilibrium potential. We assume 
it can be expressed as a function of the surface mean curvature only 




(5) 




(6) 
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where H denotes the mean curvature, / is an unknown function of H and whenever we are 
in two dimensions r = (x, y) . The presence of the square root terms in Eqs. ([s]) and ([6]) 
describes growth along the normal to the surface. 

We will further assume that function / can be expanded in a power series 

f(H)=Ko + K,H+^H' + --- , (7) 

of which only the zeroth, first, and second order terms will be of relevance at large scales. 

The result of the minimization of the potential ^ leads to the equation (see sec- 
tion III.A.3 in [26]) 

dth = /iV^/l + A [{d,,h){dyyh) - {d.,yhf] 

to leading order in the small gradient expansion, which assumes |V/i| ^ 1. Here /x = Ko, 
X = 2Ki and u = K2. We note that this equation can be expressed as the divergence of a 
current dth = V ■ j H-^'-^'^'' for j = ijVh+ | coi{D^h) ■ Vh — z/VV^/i provided some regularity 
conditions on h are assumed and where cof{D'^h) is the cofactor matrix of the Hessian 
matrix (see Eq. Q below). 

There is another way of finding the first two terms of this equation. The Hessian matrix 

(D^h) = (9) 

ydyxh dyyh J 

encodes all the information about the convexity and concavity of the interface. This matrix 
has exactly two tensorial invariants: its trace and its determinant. These are in fact the 
first two terms of our equation. A related viewpoint was adopted in [28] in the derivation 
of a model for amorphous thin film growth. 

The term proportional to the Laplacian is associated to the effect of gravity and desorption 
in the framework of nonequilibrium growth. These effects are negligible in molecular beam 
epitaxy growth, and this is precisely the context in which the VLDS equation was derived [T^ 
[16] . From a theoretical point of view, the presence of this term in Eq. ([s]) makes the 
system evolve towards the linear Edwards- Wilkinson fixed point. So theoretically it is more 
interesting to neglect this term and this way find a genuine nonlinear dynamics. Setting 
/i = in Eq. g yields 

dth = X [{dxxh){dyyh) - [dxyhf] - uV^h + e^i''). (10) 



This equation can be considered as an alternative to both SSG and VLDS equations [29] . 
Indeed, it is identical to both except for the nonlinearity, that shares the same dimensional 
analysis properties and the fact that it is conserved. 




FIG. 1: Numerical solution of Eq. (10) in the deterministic limit D = 0. The values of the 
parameters are A = 1 and ly = 0.1. We have assumed periodic boundary conditions and the 
symmetric initial condition ho = sin(27rj;) sin(27ry). One sees that holes disappear and mounds 
form as a consequence of the action of the nonlinearity. Red and orange colors denote negative 
values of h. 



In order to clarify the meaning of the nonlinear term in Eq. (10) we numerically solve 
this equation in the deterministic limit D = 0, see Fig. [Tj We see that holes disappear and 
mounds form as a consequence of the action of the quadratic nonlinearity. This phenomenol- 
ogy can be explained by means of a simple geometric argument. The surface Gaussian cur- 
vature is given in the small gradient limit by /C ^ h^^hyy — /i^^, so the dynamics favors the 
growth of those parts with a positive Gaussian curvature, which is precisely the observed 
effect. 

As a final remark, we note that the growth of perturbations with a defined wavelength 
may be achieved by means of considering the full Eq. (|8| with /i < 0. The effect of this 
type of linear instability on the VLDS equation has already been considered [30]. It could 
be possible to perform a renormalization group analysis in this case as well. One can use 
the studies on the noisy Kuramoto-Sivashinsky equation which use this technique 
as a benchmark to this purpose. 
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IV. RENORMALIZATION GROUP ANALYSIS 



Herein we will analyze the equations derived in the previous section as well as some pos- 
sible variants. Our main tool will be the dynamic renormalization group (RG) as developed 
in HIH]. 



A. Methodology 

As a first step we will explain how to implement the RG analysis in our current situation. 
Due to the particular form of the nonlinearity, the present analysis requires considering some 
subtleties which are not present in the case of more studied nonlinear terms. 

We begin with the equation 

dth = \{h^.hyy - hly) - UV'h + ar,t), (11) 

forced by a white Gaussian noise characterized by its two first moments: 

(e(r,t)) = 0, (12) 
{ar,mr',t')) = D6{r-r')S{t-t'). (13) 

Next we Fourier transform this equation to find 



h{k, ui) = hQ{k, uj) + X Go{k, uj) x 
"1 2 2 1 



xh{k-q,co-Q)h{q,Q)^, (14) 



where 



ho{k,uj) = Go{k,uj)^{k,uj), (15) 
Go{k,uj) = {-iu + ukY\ (16) 

is the Fourier transformed white noise, k = \k\ and k = {kx,ky) (the same holds for q). 
This equation is to be solved iteratively in the vertex. The resulting Feynman diagrams are 
identical to the ones represented in [5j and [H]. The dynamic renormalization group tech- 
nique |1] calculates the intermediate values of the renormalized parameters by integrating 
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out fast modes in the momentum shell e~^A < k < A, where A is the momentum cutoff and £ 
parameterizes the change across different scales. The remaining slow modes {k < A) are re- 
stored to full momentum space by means of rescaling space and time k ' = k, u' = e~^^ u, 



and field h' = e"^ h. The scaled field h' obeys Eq. (14) as well provided the renormalized 
coefficients satisfy the renormalization group flow equations, that we will calculate in the 
following. 

In the present case the noise amplitude is D, the propagator is 

1 



-iu + vk* ' 



(17) 



and the vertex is 



^ / + - k.kyq.qy^ = X j ]^{k . q^^ ^ (18) 

where q-^ = {qy, — g^)- The last equality is advantageous because it will allow us to integrate 
over the angle formed by k and q-^. If this equality were not present we would have to 
carry out a more complicated calculation reminiscent to that described in the Appendix. 
Note that here we are explicitly using we are in c? = 2. In two dimensions the situation is 
particularly simple because for any given non-zero vector there are just two vectors of a given 
modulus perpendicular to it. Our selection is arbitrary: choosing the opposite possibility 



for q does not modify the result as the scalar product on the right hand side of Eq. (18) 
is squared. Designing the optimal vertex form for a similar model built in d > 2 is a more 
technical issue. 

The first diagrammatic contribution to propagator renormalization is 



1 



87r3 {-iuj + uk^) 



■'■ \2„2 I -'- /T „ \2„2 



{k^ - qx) qy + :^{ky- qy) q^ - {k^ - qx){ky - qy)qxqy 



X 



It is important to realize that we have the simplification 
\{k^ - qxfql + \{ky~ qyfql - {k^ - q^){ky - qy)q^qy = \klql + ^fc^^^ _ k^kyq^q^. (20) 



Now using Eq. (18) we reduce this diagrammatic contribution to 

D f k^q^ cos\e) 



f kYcos\9) 



where 6 is the angle formed by k and g* , and we have carried out the hmit k,uj — t- 0. 
We note this contribution is relevant because the integrand is 0(/c^) and the propagator 
renormalizes at this same order. Diagrammatically this one loop order is identical to the 
one represented, for instance, in figure 2(a) in [Sj and figure 1(b) in |14j . 
The first diagrammatic contribution to noise renormalization is 

S^]^ j (^,j2q8^Q2y (^2^lli + 2^lll - dq dQ, (22) 

in the limit k,u ^ 0. This contribution is irrelevant because the integrand is 0{k'^). The 
corresponding diagrams at one loop are the same as those represented, for instance, in figure 
2(b) in [5] and figure 1(c) in [H]. 

The first diagrammatic contribution to vertex renormalization is irrelevant for exactly 
the same reason. Diagrammatically this one loop order is identical to the one represented, 
for instance, in figure 2(c) in [3] and figure 1(d) in [13]. In this case all diagrams have three 
external legs: one corresponds to an ongoing ki momentum and the other two to outgoing 
/c2 and ki — k2 momenta. Consider for instance the intermediate vertex corresponding to an 
ongoing leg ki and two outgoing legs q and ki — q: 

^{ki^ - qxfql + ^{kiy - qyfql - {ki^ - q^){kiy - qy)q^qy (23) 

which is of order 0{k'^). The other two intermediate vertices are trivially of order 0{k'^), 
so the order of this diagrammatic contribution is 0{k^). Automatically this implies this 
contribution is irrelevant in the hydrodynamic limit A; — )■ 0, because the vertex renormalizes 
at order O(fc^). 

The renormalization group flow equations then read 

dv 



di ^ 
dD 



3 

z -A- 



(24) 



6Atc 

D[z-2-2a\, (25) 



di 

f^=X[z + a-A]. (26) 

These equations yield the following values of the critical exponents: a = 2/3 and z = 10/3. 
The exponents are exactly the ones found for the VLDS equation at one loop order. The 
difference in this case is that the result is exact. We note that noise does not renormalize at 



any order in the loop expansion yielding the exactness of Eq. (25). This is a well known fact 



due to the conserved character of the drift of the stochastic growth equation. In the present 



case, Eq. (26) is exact too due to the fact that the vertex does not renormahze at any order 
in the loop expansion. All diagrams contributing to vertex renormalization at any order are 
irrelevant in exactly the same way we have shown herein for the one loop order diagrams. 

As we have seen, the exactness of the scaling relation a + z = 4 holds in this case. 
The current situation is genuinely different from the one concerning both SGG and VLDS 
equations. In that case vertex non-renormalization was fulfilled at one loop due to the 
cancellation of all diagrams contributing at this order. Janssen showed that the cancellation 
of diagrams contributing to vertex renormalization at two loops is not exact, what leads 
to the inexactness of the scaling relation [17] . In the present case this argument does not 
apply because there is no cancellation of diagrams: instead all diagrams contributing to 
vertex renormalization are irrelevant. So the reason underlying the exactness of the scaling 
relation is fundamentally different. 



It is important to note we have used a symmetrized vertex given by Eq. (18) after using 



the simplification Eq. (20). This choice is by no means unique, as we could have chosen a 



different vertex form, such as 



A j [{K - Qxfql - {kx - qx){ky - qy)qxqy\ , (27) 
or an infinite number of different possibilities. Among all of them, the only one which is 



symmetric under a simultaneous change of coordinates {k^ < — )■ ky.q^ < — > qy) is Eq. (18). 
Contrary to what we have shown herein, using an asymmetric form results in the dependence 
of the RG flow equations on the polar angle in Fourier space 0. Indeed, the substitution 



kx = k cos(0) and ky = k sin(0) in Eq. (27) does not lead to a vertex form independent of 
(f). Formally proceeding, this dependence is translated to the flow equations. This leads to a 
perfectly valid result whenever the resulting exponents are independent of this angle, see |29] . 
However, in the general case, some dependence on this angle may arise. The presence of this 
angle in the RG flow equations signals the generation of an integral operator in the stochastic 
growth equation when its coefficients become renormalized. This consequently leads to a 
more complicated mathematical treatment, as it implies the necessity of restarting the RG 
procedure including the newly generated operators. In order to avoid this, the symmetric 
form Eq. (18) must be used. We discuss in Appendix |A] the different results which may arise 



from the RG analysis employing asymmetric vertex forms. 
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B. Results for different models 



Now we will explore different stochastic growth equations with the methodology explained 
in the previous section. Our goal is getting a deeper understanding of the properties of the 
nonlinearity we have considered so far or variants of it. In particular, we will show that vertex 
non-renormalization is present in those cases in which noise renormalization is present. 



We start with a variation of Eq. ( 11 ) 



dth = \{Kxhyy - hi ) - vV^h + ^(f, t) 



(28) 



where now the noise is conserved, i. e. it is white and Gaussian and its two first moments 
are given by 



W,t)) 

mmf\t')) 



0, 



-DW'^5{f-f')5{t-t') 



(29) 
(30) 



As we have already outlined, this equation can be considered as an analogue of the SSG one. 
The renormalization group fiow equations can be obtained following analogous arguments 
to those of the previous section and read 

3 \^D' 



dv 

M 

dD 

~d£ 
dX 



z-A- 



6Att 

D[z- A- 2a], 
\[z + a-A]. 



(31) 

(32) 
(33) 



From here we get the critical exponents a = and z = 4, this is, the SSG universality class. 
We emphasize that the last two RG fiow equations are exact, i. e. both noise and vertex do 
not renormalize at any order in the loop expansion, for exactly the same reasons as in the 
previous section. 

As the next example we consider the following equation 



dth = X{hxxhyy - hly) + vV^h + i{r, t) 



(34) 



where the noise is again white and Gaussian and its first and second moments are given 



by Eqs. (29) and (30) respectively. We note that the sixth order differential operator has 
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already been deduced in the context of nonequilibrium growth from atomistic models [33] . 
The renormalization group flow equations become 



di 

dD 

~dl 
dX 



z — 6 



45 X^D 



1287r 1/3 
D[z-4:-2a], 

X[z + a-A\. 



We find again the same exponents a = and z 
renormalization (both results are exact again). 



(35) 

(36) 
(37) 

4 due to noise and vertex non- 



The next equation that deserves our attention is a different variation of Eq. (11) 

/iL) -z/V'/i + e(r,t), 



dth = XiK^hyy - hiy) - + e(r, t), (38) 

where the noise is as always white and Gaussian but this time its mean and correlation are 

(e(r,t)) = 0, (39) 
mmr'.t')) = DV%f- r')5{t - t'). (40) 

This sort of higher order contribution to conserved noise has been previously considered in 
the literature [31]. In this case, a simple dimensional analysis shows that the nonlinearity 
is not relevant and leads to the exponents a = —1 and z = A. In this context a negative 
exponent a is identified with a fiat interface. 



In the case of Eq. ( 34 ) 



dth = X{h^,hyy - hly) + z/V^/i + i{f, t) 



(41) 



where the noise is white and Gaussian and its mean and correlation are given by Eqs. (39) 



and (40) respectively, the RG analysis yields a negative exponent a as well. This can be 



read from the following RG flow equations: 

dy 
di 
dD 



V 



Z-6 + 



45 X'^D 
1287r^ 



D 



z-Q-2a + 



3 X^D 
5127r 1/3 



dX 
d£ 



X[z + a-4] 



(42) 
(43) 
(44) 



Thus the exponents become a = 240/179 - 2 ^ -0.66 and z = 6 - 240/179 ^ 4.66. Note 
that vertex non-renormalization is again present and exact, while the noise is renormalized 
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at one loop order. This result, noise renormalization in a conserved model, was already 
found in the literature for this kind of noise and a different drift 
Now we direct our analysis to higher order equations like 



dth = \{K^hyy - hi ) - vV^h + e(r, t) 



(45) 



where the noise is white and Gaussian with mean and correlation given by Eqs. (29 ) and (30 ) 



The eighth order differential operator can be considered as a higher order diffusion mech- 
anism in the context of nonequilibrium growth, when we regard continuum equations as 
hydrodynamic descriptions of atomistic models [33]. This equation can be analyzed with 
the RG technique and we have found that neither noise nor vertex do renormalize at any 
order, what yields the exact exponents a = and z = A. 
The previous equation 



dth = \{hxxhyy - hi ) - vV^h + e(r, t) 



(46) 



can be studied with a Gaussian white noise having as first moments Eqs. (39) and (40). Its 



analysis by means of the RG method yields a negative value of exponent a. In this case the 
flow equations read 

' . 1 X^D' 



dv 
dl 
dD 
~d£ 
dX 
di 



D 



z-Q-2a + 



3 X'D 
128n 1/3 



X[z + a-A] 



yielding the exponents z = 1768/381 ~ 4.64 and a = 4 



(47) 
(48) 
(49) 

-0.64. Again we find that 



vertex non-renormalization is exact and noise renormalizes at one loop order. 
Our last example is the following equation 



dth = XV^{hxxh 



yy 



Ky) 



uV^h + ^{f,t), 



(50) 



where the white and Gaussian noise is characterized by the moments Eq. (39) and (40). 



We found this stochastic differential equation specially interesting because its RG analysis 
is meaningless as neither propagator, nor noise, nor vertex renormalize at any order in this 
case. So a different technique should be employed for this kind of model. One possibility 
would be the use of self-consistent expansions [T81423] as we mentioned in the Introduction, 
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but we have not explored this possibihty yet. We also note that, due to the high order of 
differentiation in every term of this equation, its numerical analysis should be a hard task 
too. 



V. ORIGINS OF SCALING 

We now examine the physical origin of scaling relations. The exactness of the relation 
a + z = 4 in all our equations can be traced back to expansion ([t]). Higher order terms are 
meant to describe smaller scale physics. Accordingly, the exactness of this scaling relation 
indicates that the nonlinearity totally dominates over surface diffusion in the large scale. 
Considering the full Eq. ^ with /i > we find, in agreement with Eq. ([T]), that the 
Laplacian dominates. So expansion ([t]) provides a simple geometric interpretation of the 
renormalization group results and the exactness of this scaling relation. 

The exactness of the exponents for the SSG and VLDS universality classes, found in our 
equations, was due to noise and vertex non-renormalization. Noise non-renormalization is a 
common fact, and its origin is in the conserved character of the drift, which has been shown 
in previous works [HI [161 [H] • The case of the vertex is new. Indeed, the vertex structure 

X[l(k-q^y (51) 



2 

implies that all contributions to the renormalization of the vertex are 0{k^), while the ver- 



tex renormalizes at 0{k ). The concrete calculations have been carried out in Sec. IV A 
In other words, the vertex is not renormalized because all the generated Feynman diagrams 
correspond to irrelevant contributions. This is easily seen by noting that in the renormaliza- 



tion group calculation the vertex contributions result from integrating against measure (51). 
The procedure is identical to the one carried out for the KPZ equation [5J mutatis mutandis. 
Any possible contributions to the renormalized coefficients come from integrands that are 
exactly 0{k^), as higher order contributions are irrelevant. It is immediate realizing that 



the 0{k^) contribution of measure (51) is identically zero. In consequence, noise and vertex 
do not renormalize at any order in the loop expansion. This is in contrast to what happened 
to the SGG and VLDS equations, for which the vertex one loop diagrams were relevant but 
they cancelled each other out [HI [16] . 

While the renormalization group is a perturbative technique which may allow access to 
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the critical exponents, it is also possible to unveil (at least part of) the critical behavior by 
means of finding the stationary probability distribution of a given Langevin equation. This 
approach has being successfully exploited in the past mainly (but not only [37j) in 

those cases in which the stationary state is Gaussian. In our case it is possible as well, as 
was shown in [29], to formally find the exact stationary probability density corresponding 



to Eq. (|ll|) using its gradient flow structure 

V[h] = 



6V 

dth = -—+ar,t), 



XKhyh^y - -{y'^hf 



dx dy. 



(52) 



Then the probability density reads 

V ~ exp 



XKhyh^y - -{y'^hf 



dx dy 



(53) 



and it is evidently not Gaussian. This distribution is of course not normalizable due to 
the presence of the cubic term, and consequently this expression is purely formal. Its form 
suggests that the interface profile could be initially dominated by the quadratic term, and 
after some transient the noise could drive the interface out of this potential well. Then 
for long times the cubic term would dominate, as the dynamics would be separated from 
the initial Gaussian behavior by the nonlinearity. This simple picture, derived from the 
formal stationary probability distribution, agrees with the results of the renormalization 
group analysis. 

VI. CONCLUSIONS 



In this work we have shown that the scaling relation a + 2; = 4, which originates in the 
vertex non-renormalization and which is inexact in the renormalization group sense for both 



SSG and VLDS equations, becomes exact in the case of Eq. (10) with both conserved and 



non-conserved noises. The physical origin of this scaling relation is not the existence of 
a symmetry (such as a functional Galilean invariance) but the existence of the geometric 
principle described by Eqs. (|5|-([6])-([7]). 

We have proven the exactness of this scaling relation in = 2, in the renormalization 



group sense, for a series of stochastic differential equations presented in Sec. IV B Therein 



we have found that vertex non-renormalization is exact for the considered nonlinearity in- 
dependently of whether noise is renormalized or not. The validity of these results extends 
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to all examined models, except the last one. In that case, the higher order of differentia- 
tion precludes the access to the critical behavior by means of the employment of the RG 
technique, as we have shown. 

We have focused on conserved surface growth but the existence of scaling laws happens 
as well in non-conserved cases. Notably the KPZ equation obeys the relation a + 2 = 
2. Although this scaling relation has been traditionally attributed to the existence of a 
Galilean invariance in the equation, recent numerical work suggests that it is possible to 
find discrete schemes which do not obey this symmetry and are still able to reproduce the 
KPZ universality class [SHU]. Perhaps the existence of a geometric principle, akin to the one 
discussed here, could be a possible alternative explanation. The origin of such a principle is 
not clear so far, but the variational formulation of the KPZ equation [38] might have some 
relation with it, in case it exists. 
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Appendix A: RG Analysis with an asymmetric vertex form 



As mentioned in section IV A using asymmetric vertex forms leads to the explicit appear- 
ance of the polar angle in Fourier space. This fact signals the generation of new operators 
and requires the employment of a more sophisticated mathematical treatment that will be 
outlined at the end of this appendix. The aim of this appendix is showing that, nevertheless, 
it is possible to formally proceed with the calculations ignoring the presence of this angle 
and still find the good result. Sometimes this is possible, sometimes it is not. Herein we 
exemplify these different possibilities. 

We consider the equation 

dth = X{K,hyy - hly) - UV^h + (Al) 

First of all we shall concentrate on the non-conserved noise case ^ = We employ the 
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following asymmetric vertex form 

h{k, u) = ho{k, u) + \ Go{k, u) x 

[{ka, - Qxfql - {kx - qx){ky - Qy) Qx Qy] 

dfl dq 



where 



X h{k - q,uj -Q) h{q, Q.) 

ho{k,iu) = Go{k,iu) ^{k,iu), 
Goik,uj) = {-iu + uk'^)-\ 



(27r)3 ' 



(A2) 

(A3) 
(A4) 



^ is the Fourier transformed white noise, k = \k\ and k = {kx,ky) (the same holds for q). 
The renormalized coefficients satisfy the RG flow equations, which at one loop order are 
given by 



di 
dD 

H 

dX 



— V 



4 - - — ^cos (0) 



87ri/3 
D[z-2-2a], 

A + a - 4] . 



(A5) 

(A6) 
(A7) 



The variable is the polar angle in Fourier space. An alternative way of thinking about 
this angle is through the relation 



cos (0) 



(A8) 



^xx ~t~ ^yy 

where denotes the spatial Fourier transform. The values of the exponents are independent 
of this angle and take the correct values a = 2/3 and z = 10/3. 
Now we move to the case of conserved noise, i. e. to the model 



dth = X{hxxhyy - hi ) - + ^ 



(1) 



(A9) 



Repeating the dynamic renormalization group calculations we arrive at the flow equations 
for the effective parameters 



dp 
dl 
dD 

~dl 
d\ 



D 



87ri/3 



^ - 4 - 2q; + 



COS^(0) 

647rz/3 



X[z + a-4\ 



(AlO) 
(All) 
(A12) 
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The resulting exponents are independent of and consequently they are the correct ones: 
a = and z = 4. Note however that the noise renormalization that arises in this case should 
be interpreted as a spurious result since it does not appear in the case of a symmetric vertex 



form, cf. Eqs (31 )-(32)-(33). 



Now we move to the last model 



The renormalization group leads us this time to the flow equations 
dv 

di ^ 



(A13) 



dD 
~di 
d\ 
M 



V 



D 



6 + 



cos2(0){19- lOcos(20)} 



6477 z/^ 



So we find the 



X[z + a-A]. 
-dependent exponents 



(A14) 

(A15) 
(A16) 



8 



a 



-4 + 3cos2(0){19- 10cos(2(/))}' 
8 

z = A + 



(A17) 
(A18) 



4-3cos2(0){19- lOcos(20)}' 

In this case these values have to be considered invalid because it is not possible to recover 
the correct exponents found in the analysis performed with the symmetric vertex form from 



Eqs. (A17, A18). Indeed, averaging out from either Eqs. (A14, A15, A16) or (A17, A18) 



does not reproduce the actual values of the exponents a = and z = A known form Eqs. (35 



36, 37). The correct way to proceed now would be deriving the stochastic growth equation 



with renormalized coefficients. This leads to the appearance of new operators in the equation. 
An example with a fourth order equation can be found in [39]. In this example, contrary to 
ours, a new term in the equation is generated; in our present case the differential operators 
become integro-differential. In fact, cos2(0) gives rise to a linear integral operator which in 
Fourier space reduces itself to a multiplication by the factor 

(A19) 



1.2 I 1.2 
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of the field over which it is apphed |10] . The resulting stochastic growth equation should be 
analyzed again by means of the RG technique. Obviously, this is much harder than using a 
symmetric vertex form. 

In summary, it is possible to analyze stochastic growth equations using asymmetric vertex 
forms. One can formally proceed, ignore the presence of polar angles in Fourier space, and 
end up with 0— dependent flow equations. If the resulting exponents are independent of 
these angles one can conclude the calculation. If there is (f) dependence in the exponents, 
then one has to derive the equation with the newly generated operators and restart the 
analysis. The overall procedure is more complicated than using a symmetric vertex form. 
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